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1. INTRODUCTION 

Currently, researchers pay much attention to 
graphene-based superlattices. In particular, methods 
of molecular dynamics were used to calculate graphene- 
based superlattices with periodically located rows of 
vacancies [1]. Then superlattices of single-atom thick¬ 
ness, formed by lines of pairs of hydrogen atoms ad¬ 
sorbed on graphene, were calculated within the den¬ 
sity functional theory [2] . Rippled graphene, which can 
be considered as a superlattice with a one-dimensional 
periodic potential of ripples, was investigated in [3, 4]. 
Analytical studies were performed on superlattices ob¬ 
tained by applying a periodic electrostatic potential [5- 
8] or periodically arranged magnetic barriers [9-11] to 
graphene. A detailed review of graphene-based semi¬ 
conductor heterostructures can be found in [12]. 

We investigate planar superlattices based on gap¬ 
less graphene and its gapped modifications. The main 
concepts of these superlattices were formulated in [13], 
where a dispersion relation for charge carriers in these 
structures was derived. Then a very simple example of 
such superlattices was considered: two-type superlat¬ 
tices composed of alternating strips of gapless graphene 
and its band gap modification. 

In [14] we studied graphene-based quantum wells 
(QWs): planar heterostructures based on graphene, 
where gapped graphene modifications play a role of 
potential barriers. In particular, it was shown that 
the energy spectrum of asymmetric QWs (containing 
different gapped graphene modifications) is split with 
respect to pseudospin: the dispersion curves in dif¬ 
ferent valleys do not coincide. This result suggests 
that pseudospin splitting should also occur in polytype 
graphene-based superlattices. This splitting is similar 


in many respects to the spin splitting of the energy 
spectrum in narrow-gap heterostructures [15, 16]. 

Note also that recently we proposed and investi¬ 
gated another version of planar superlattices, which is 
peculiar for graphene [17]. Using the Fermi velocity en¬ 
gineering in gapless graphene, one can fabricate struc¬ 
tures with periodically modulated Fermi velocity. 

In this paper, we discuss the pseudospin splitting 
of the energy spectrum of a polytype superlattice. 
As an example of polytype graphene-based superlat¬ 
tices, we consider a three-type superlattice in the form 
A — B — C, where A and C are gapped graphene modi¬ 
fications with different bandgap widths, and B is gap¬ 
less graphene. This is the simplest example of a poly¬ 
type superlattice which retains fundamental features of 
polytype superlattices. 

To implement a three-type graphene-based super¬ 
lattice, we purpose a version involving gapless graphene 
and its gapped modifications obtained as a result of 
deposition of a gapless graphene sheet on a particular 
substrate (e.g., hexagonal boron nitride h-BN) or due 
to deposition of particular atoms or molecules (for ex¬ 
ample, CrOa) on the surface of gapless graphene. This 
version is presented pictorially in Fig. 1. 



Fig. 1. A version of a three-type superlattice: a graphene 
sheet on a stripped substrate consisted of SiOa with de¬ 
posited stripes of CrOa molecules and h-BN. 
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2. MODEL 


Let region I be a layer of gapped graphene modifi¬ 
cation with a bandgap half-width A/ and work func¬ 
tion V/; the thickness of this region is dj. Region II is a 
layer of gapless graphene with a thickness djj and zero 
work function. Region III is a layer of gapped graphene 
modification with a bandgap half-width A///, work 
function V/ 77 , and thickness dm (Fig. 2). The su¬ 
perlattice period is d = ^7 -I- dn -I- dm. 

Expecting to find pseudospin splitting in the ener¬ 
gy spectrum of the superlattice under consideration, 
we should investigate the solutions to the 4x4 ma¬ 
trix Dirac equation describing both valleys. As was 
shown in [18], one can pass to the 2 x 2 matrix equation 
for the components of the eigenspinor of pseudo-parity 
ope-rator. Note that this equation contains explicitly 
the operator eigenvalue — pseudo-parity A, which dis¬ 
tinguishes states from different valleys: A = -1-1 for 
the states near the K point and A = — 1 for the states 
near the K' point. We will use specifically this equa¬ 
tion to make the computation less cumbersome: 


{vfctxPx + XvFCTyky + Ao-j -I- V) ip\{x) = Extpxix), 

( 1 ) 

where vf is the Fermi velocity (it is constant in all 
superlattice regions); <Jy, and are the Pauli ma¬ 
trices, which act in the sublattice space; Px = —idx 
is the X component of the momentum operator (from 
here on, Planck’s constant is assumed to be unity); 
ky is the y component of the quasimomentum (charge 
carriers move freely along the y axis); and A and V are 
quantities periodically changing along the x axis with 
a period d, 



Fig. 2. Energy band diagram of a three-type superlattice, 
supercell of which consists of a gapless graphene strip and 
two strips of its gapped modifications. The position of one 
the lower minibands of interest is shown by a gray strip. 


A = 


V = 


{ Aj, (n — l)d < X < (n — l)d -I- dj, 

0, (n — l)d + di < X < nd — dm, 

A 777 , nd — dm < x < nd; 

{ Vi, (n — l)d < x < {n — I)d-|- dj, 

0 , (n — l)d + di < x < nd — dm, 

Vm, nd — dm < x < nd. 


where n is an integer enumerating supercells (see 
Fig. 2). 

Let us write solutions to Eq. (1) in each region 
within the nth supercell: 


(i) in region I (0 < a: < ^ 7 ), 

= ( 2 ) 

where 

^ exp(-A:ia:CT,), 

(±) _ ivF(ki ± Xky) 

“ Ex + Ai - Vi’ 


(ii) in region II {di < x < di + dn). 


II 

< 1 

f “in] 

Uv ■ 

(3) 

where 





1 

exp(ik2xaz). 


o(±) _ '^^’(^2 ± iXky) ^ 




(hi) in region III {di + dn < x < d) 

= (“PP (4) 

where 

^ -p-p ^M-ksxcr,), 

(±) _ ivFiks ± Xky) 

“ Ex + Aiii-Vm' 

A is the normalization factor for all three regions. 
The ki, k 2 , and k^ values are related by the correspon¬ 
ding expressions to energy Ex'. 

= ^Fjii ~ — Vijii)^ + v%ky, 

vlkl =El- vlkl 

The ^’^‘4 ^^Xn values are constants, which can be 
found from the boundary conditions. 

Now, based on solutions (2)-(4), we can construct 
the wave function of charge carriers in the superlattice. 
In this study, we are not interested in the interfacial 
states (see, for example, [18]); therefore, vf is assumed 
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to be the same in all regions. In this case, the boundary 
conditions for the superlattice include the requirement 
for the wave function continuity within each supercell 
and between neighboring supercells [14, 19]: 

(a) continuity of the wave function between regions 

I and II at the boundary x = dj 

i ) ( 6 ) 

\^\n / \^An / 

(b) continuity of the wave function between regions 

II and III at the boundary x = dn 

+ dii) = n<-Xidi + du) , (7) 

(c) continuity of the wave function at the boundary 
between region III nth supercell and the region I of the 
(n + l)th supercell 

=^A/ci(o) ■ ( 8 ) 

\^\n J V'^An+l/ 

In view of the periodicity of the system, the Bloch 
conditions 


3. TRANSFER MATRIX 
AND DISPERSION RELATION 


Furthermore, we will determine the transfer matrix 
(T matrix) that links coefficients of the wave function 
in neighboring supercells [13, 17]: 


0 ) \ 
An+1 I 
0 ) 

An+1 / 




( 10 ) 


Definition (10) shows that there may be three ver¬ 
sions of the T matrix linking the corresponding pairs of 
constants from solutions (2)-(4) for regions I, II, or III. 
This arbitrariness in definition is allowable, because all 
three versions are interrelated by cyclic permutation of 
D matrices, and the dispersion relation includes the T- 
matrix trace. In particular, using equalities ( 6 )-( 8 ), 
one can easily derive the expression 

^i%id) {il[%idj + dn)y' 

X {dj + du) (ffit ^di)) (di). 

( 11 ) 


The dispersion relation has the form [13, 17] 


'^xl ix + d) = (x) exp(f fc,rd) (9) 

must also be satisfied for all three regions (j = 1, 2, 3). 


Trrj^^ = 2cos(A:^d). (12) 

Having calculated the T-matrix trace, we obtain the 
dispersion relation in the form 


{A‘ 


-92 


cos(fc2d//) -k sin(fc2d//) expi-ksdin) 

cos(fc2(i//) -k sm{k2du) exp(A:3d///)| exp(-fcid/) 

g[~'^ cos(fc 2 dii) - fi ^ sin(A :2 dii) \ expiksdiii) 

'' gi~'' cos(fc 2 d//) - / 2 ~^ sin(fc 2 d//) exp(-A:3d/77)| exp(fci(i/) = 
8vpkik2k3 


+ |9( 


(+) 


Ex (Ex + Aj — Vi) {Ex + Am — Vm) 


cos{kxd), 


(13) 


where the following designations are introduced: 


9 ^^ 


9 P^ 


vpiki ± Xky) vpiks T Xky) 


■ + 


Ex + Ai — Vi Ex + Am — Vm 
VFk2 


Ex 


9i j 




VFk2 (±) 



VF{ki ± Xky) 
Ex + Ai- Vi 


VFjka ± Xky) 

Ex + Am — Vm ’ 


f(±) 

J 2 


1 \k — ( (^1 =*= _ ’’Jpiks T Xky) \ 

^ ^Ex \Ex + Ai — Vi Ex + Aui — Vm) 

. ., XF / VFjki ± Xky) VFjkz ± Xky) \ 

^ ^Ex\Ex + Ai — Vi Ex + A777 — Vm ) 


+ (fci ± Xky){k3 =F Xky) 

{Ex + Ai — Vi){Ex + Am — Vm) ’ 

+ (fci ± Xky){k3 ± Xky) 

{Ex + Ai — Vi){Ex + A 777 — Vm) 


An analysis of Eq. (13) shows that the necessary condition for pseudospin splitting of the energy spectrum 
is the absence of electron-hole symmetry. The symmetry of the system in the energy space with respect the 
replacement E —>■ —E should be absent, and V 7 or V 777 should differ from zero. A sufficient condition for 
pseudospin splitting of the energy spectrum of the superlattice is the inequality 


Af A^ 
Vi ^ Vm ■ 


(14) 
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Dispersion relation (13) can be simplified in the 
case where QWs exist for only electrons or for only 
holes [13]. In the case of electrons, one must assume 
that Vi = A/ and Vm = A///; then relation (13) 
yields 


cosh(fcid7) cos(fc2(i//) cosh(fc3C?///) 

+ i sinh(/cid/) cos{k2dii) smh{k3diii) 

+ x [~2 sinh(fcid/) sin(fc 2 d//) cosh(fc 3 (i///) 


-\-X- 


3 2 ^ cosh(fci(i/) sin(fc 2 (i//) sinh(fc 3 d///) = cos{kxd), 

(15) 


where = Xij ± x~j and Xi^j = ki/kj. After 

the analytical continuation ki —>■ iki and k^ —> ik^, 
Eq. (15) coincides with the well-known nonrelativistic 
equation (see, for example, monograph [20]). Inequa¬ 
lity (14) is not fulfilled, and pseudospin splitting of 
the energy spectrum is absent. 

The limiting transition from superlattice to indi¬ 
vidual QW occurs if the potential barriers are suffi¬ 
ciently wide. Under these conditions, kidi and k^dm 
are much larger than unity, and the dispersion relation 
for superlattice (13) yields 


tan(fc2d//) 


9i 


(-) 



(16) 


for an individual QW, the pseudospin splitting in which 
was considered in [14]. 

4. NUMERICAL CALCULATIONS 


We performed calculations for the aforementioned 
gapped graphene modifications: graphene deposited on 
h-BN strips (A/ = 26.5 meV) and graphene with de¬ 
posited Cr 03 molecules (A/// = 60 meV). The strip 
widths were taken to be di = dm = 21.3 nm and 
dji = 42.6 nm. 

For simplicity, we assume that U/ = 0. The 
exact value of the work function for the second gapped 
graphene modification is unknown; at the same time, 
according to the energy band calculations, it is nega¬ 
tive [21]. We chose Vm to be -10 meV. 

The calculation results for the lower electron mini¬ 
band are presented in Fig. 3. One can see that the 
electron valleys for different pseudospins (A = ±1) un¬ 



derwent splitting. The extremum of valley is at the 
point \ky. For the superlattice under consideration, 
we have the following parameters: ky « 1.9x 10^ cm“^, 
E\{\ky) « 13.5 meV, and energy splitting Ap^ = 
~ meV for the lower elec¬ 
tron miniband and ky « 1.7 x 10"^ cm“^, E^{Xky) « 

-15.5 meV, and A^J, = Et+i(fcJ) - i^L-i(^A) « 
0.05 meV for the upper hole miniband. The thus 
formed energy gap is Eq = E^{ky)—E^{ky) « 29 meV. 

5. CONCLUSIONS 

We proposed a new class of systems based on 
graphene: planar polytype superlattices composed of 
gapless graphene and its different gapped modifica¬ 
tions. These superlattices, formed by multiply re¬ 
peated asymmetric QWs, are characterized by pseu¬ 
dospin splitting of the energy spectrum. We investi¬ 
gated the conditions for occurrence of this splitting. 

The spectrum calculated by us is similar to that 
considered by Bychkov and Rashba [22]. The rela- 
tioship between our Hamiltonian and the Bychkov- 
Rashba Hamiltonian was considered by us in de¬ 
tail in [14]. 

The results obtained barely change if we consider 
smooth heterojunctions instead of sharp ones [14]. 

Due to the pseudospin splitting of the energy spec¬ 
trum, the bandgap in this superlattice becomes indirect 
(in contrast to the superlattice considered in [13]). The 
superlattice analyzed by us is similar to a narrow-gap 
indirect semiconductor and can be used in IR optics. 
The characteristic bandgap width is Eq = 30 meV. 
The spacing between the extrema of electron and hole 
valleys with identical X is ky = ky + ky = 4x 10^ cm“^. 
The desired parameters of superlattice spectrum (Eq 
and ky) can be obtained by changing the superlattice 
characteristics {dijijn, Vjjn, and Aijn). 
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